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Detecting flux creep in superconducting YBa2Cu307 5 thin films via damping of the 

osciUations of a levitating permanent magnet. 
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The damping of the osciUations of a small permanent magnet (spherical shape, radius 0.1 mm) 
levitating between two parallel epitaxial YBCO films is measured as a function of oscillation am- 
plitude and temperature. At small amplitudes the dissipation is found to be orders of magnitude 
lower than in bulk YBCO, Q-factors exceeding one million at low temperatures. With increasing 
amplitude the dissipation becomes exponentially large, exceeding the bulk values at large drives. 
We describe our results by calculating the ac shielding currents flowing through trapped flux whose 
motion gives rise to electric flelds. We flnd dissipation to originate from different mechanisms of 
flux dynamics. 



A permanent magnet levitating above, or suspending 
below, a superconducting surface has an equilibrium po- 
sition which is determined by magnetic forces due to 
screening currents and trapped flux [Q. Except for a 
small logarithmic long-term relaxation of the magnetiza- 
tion of the superconductor this position is stable. The 
system is in a metastable state with stationary flux lines 
and consequently without any decay of the dc supercur- 
rents, i.e., without any energy dissipation. However, os- 
cillations of the levitating magnet about its equilibrium 
position give rise to ac magnetic flelds at the surface of 
the superconductor which modulate the shielding cur- 
rents. The initial state is now disturbed. The super- 
conductor periodically tries to relax to more favorable 
metastable states by moving flux lines back and forth 
as a result of the oscillating component of the Lorentz 
forces. This leads to energy dissipation which we can de- 
termine from measurements of the oscillation amplitude 
of the magnet as a function of an external driving force. 

In our present work we have investigated the energy 
dissipation when the magnet is levitating between two 
horizontal YBa2Cu307_5 (YBCO) epitaxial thin fllms. 
We flnd that at small amplitudes and low temperatures 
the damping is several orders of magnitude smaller than 
with YBCO bulk samples studied in earlier work 
leading to Q-factors above 10^. This result may be of 
significance when application of superconducting levita- 
tion, e.g., for micromechanical bearings, is considered [Q. 
At large driving forces and at temperatures above 77 K, 
however, the dissipation grows exponentially with the os- 
cillation amplitude and ultimately even exceeds the bulk 
values. We discuss our results in terms of thermally acti- 
vated flux flow phenomena in the superconducting fllms. 

Our experimental method has been described in de- 
tail in earlier work on sintered YBCO [^, or niobium 
1^, and also when using the oscillating spherical mag- 
net for hydrodynamic experiments in superfiuid helium 
W. It consists of placing a magnetic microsphere made 



of SmCo5 (radius 0.1 mm) inside of a parallel plate ca- 
pacitor (spacing d — I mm, diameter 4 mm) having elec- 
trodes which in our present work are made of YBCO 
epitaxial thin films (thickness 6 = 450 nm and in a sec- 
ond experiment 190 nm) laser deposited on insulating 
substrates ||^. We prepared the YBCO films on SrTiOs 
(lower electrode) and Y203-stabihzed Zr02 (upper elec- 
trode) substrates cut with an inclination of 2 degrees off 
the (001) orientation. The inclination improves both flux 
pinning and, due to initiating a terrace growth, crys- 
tallinity of the films. The lower electrode is protected by 
a 20 nm thick PrBa2Cu307 epitaxial layer. The YBCO 
fllms have an extremely sharp superconducting transition 
at Tc = 89.5 K measured by an ac susceptibility method, 
the complete transition width being 0.1 K. Before the 
capacitor is cooled through Tc we apply a dc voltage of 
about 800 V to the bottom electrode. Therefore, the 
magnet carries an electric charge q of about 2 pC when 
levitating. The dc voltage is then switched off and oscil- 
lations of the magnet can be excited with an ac voltage 
Uac ranging typically from 0.1 mV to 20 V and having 
a frequency at the resonance of the oscillations (« 300 
Hz). These vertical oscillations induce a current qv/d 
in the electrodes, where v is the velocity of the mag- 
net, which is measured by an electrometer connected to 
a lock-in amplifler. For a given driving force F — qUac/d 
we measure the maximum signal when slowly sweeping 
the frequency towards resonance. We then vary the driv- 
ing force at constant temperature. In order to infer the 
velocity amplitude from the measured current amplitude 
and the driving force from the ac voltage we determine 
the charge q by recording a resonance curve and analyz- 
ing its width (which is given by the damping coefficient) 
and its height which then determines the charge. The 
driving forces range from 10~^^ N to 10~^ N, the ve- 
locity amplitudes from 0.1 mm/s to 50 mm/s, and the 
oscillation amplitudes a = v/uj from 50 nm to 30 fim. 
The dissipated power varies from 10^^ W down to below 
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10~^ W which is less than what is usually dissipated 
when investigating flux creep by measuring the ac sus- 
ceptibility or the current- voltage characteristic of a type 
II superconductor. 




driving force amplitude [nN] 
FIG. 1. Resonant velocity amplitude of the oscillating 
magnet as a function of the driving force at various temper- 
atures. Note the steep increase at small driving forces and 
a logarithmic dependence at large drives. The lines are fits 
of our model to the data, see text. A typical linear depen- 
dence observed with melt-textured bulk samples is shown for 
comparison. 
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FIG. 2. Velocity amplitude at very small driving forces in- 
dicating the linear initial increase at various temperatures. 
The lines are fits of our model to the data, see text. 

A summary of our data obtained with the 450 nm films 
above 77 K is shown in Fig. 1. At small driving forces 
we find a steep initial increase of the velocity amplitude 
of the magnet, whereas at driving forces above ca. 10 nN 
there is only a logarithmic dependence. This implies an 
exponential increase of the dissipation. At lower tem- 
peratures this logarithmic regime could not be reached 
because the oscillations became unstable at large ampli- 



tudes. The small-signal behavior is shown in Fig. 2 with 
a largely expanded abscissa. Only at very small ampli- 
tudes we do find a linear regime between the driving force 
amplitude F and the velocity amplitude u, i.e., F = jv 
with a strongly temperature dependent coefficient 7(T), 
see Fig. 3. While 7 decreases considerably below the crit- 
ical temperature Tc, it obviously levels off at 4 K. Here 
the Q-factor ma;/7 (m = 5.2 • 10~® kg is the mass of the 
magnet) exceeds 10^ The steep drop below Tc follows 
approximately a [Tc — T)~^ law as can be seen in the in- 
sert of Fig. 3. Comparing these results with our data on 
bulk YBCO 1^ we note that at 78 K 7 has decreased by 
three orders of magnitude. 

For a more quantitative analysis we first calculate the 
sheet current on the surface of the superconductor by as- 
suming that the magnetic field on either surface is given 
by the dipolar fields of the magnet and the first image 
dipole, neglecting images further away and also modifi- 
cations of the field due to trapped flux lines. From the 
remanencc of the magnet and its radius we calculate its 
magnetic moment to be 4 • 10^® Am^. We further as- 
sume that the magnet levitates in the middle between 
the electrodes, i.e. at a height h = ho = 0.5 mm 0. 
It is simple to calculate the dc sheet current on the sur- 
face for a particular orientation of the dipole with re- 
spect to the surface The orientation of the dipole 
parallel to the surface has the lowest energy but because 
of trapped flux we cannot rule out that some other ori- 
entation applies. Although the variation of the surface 
current J — {Jx{x,y,hQ),Jy{x,y,hQ)) [x and y are co- 
ordinates on the surface) is quite complex and changes 
drastically with the orientation of the dipole, the maxi- 
mum value Jmax is fairly insensitive to the orientation: 
Jmax = 4400 A/m when the dipole is perpendicular and 
Jmax = 5100 A/m for the parallel orientation. In either 
case the magnetic field is below i?ci, but some flux re- 
mains trapped because we have a field-cooled situation. 
We then calculate the ac sheet current AJ due to the 
oscillations of the magnet about its position at ho: 

where a is the oscillation amplitude. Again the current 
distribution is complicated but the maximum amplitude 
of A J is insensitive to the orientation of the dipole. An 
oscillation amplitude a = 0.5 fim corresponding to a ve- 
locity amplitude of f « 1 mm/s leads to a maximum 
current amplitude of 14 A/m for the perpendicular ori- 
entation and of 15 A/m for the parallel one. The current 
is localized in a small area of the superconducting sur- 
faces below and above the sphere. This ac sheet current 
is actually an oscillating current density j in the film 
which exerts an oscillating Lorentz force on the trapped 
flux lines. Some flux lines will be set into motion thereby 
creating an electric field E = Vfn^o, where Vf is the flux 
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line velocity and n is the number of flux lines per area. 
Because j and E are parallel the dissipation per unit vol- 
ume is given hy j ■ E. If, for simplicity, we use an average 
current density (j) inside the film of thickness (5, we can 
describe the dissipation per unit area by AJ ■ E, where E 
depends on (j) = AJ/6 We can calculate the station- 
ary oscillation amplitude at resonance by employing the 
energy balance between the gain from the driving force 
and the loss per cycle (period t): 



F-vdt 



AJ ■ E dxdydt 



(2) 



J x^y 



In the linear regime F = -fv the loss (r.h.s. of Eq. H) can 
be described by a surface resistance Rg oc 7, namely 



Rs ■ (A J)2 dxdydt = R, 



x,y 



[Ajf dxdydt. 



^ x^y 



(3) 



In Fig. 3 the Rs values are given for the parallel orien- 
tation of the dipole . From the surface resistance we 
can infer the imaginary part of the penetration depth of 
the ac field being given by A" = Rs/fiQUj. In our case 
(w/27r « 300 Hz) we find A" to drop from 84 nm at 85 K 
down to only 0.2 nm at 4 K. While at present we do not 
have a model for a quantitative analysis of these data, 
we note that a steep decrease of the surface resistance 
with temperature is also observed at both radio and mi- 
crowave frequencies and is believed to be due to small 
oscillations of the trapped flux lines or, in zero field, due 
to normal fluid electrons. 
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FIG. 3. Linear coefficient 7 = F/v from Fig. 2. Insert: 
same data (ordinate unchanged) plotted vs Tc — T. The sur- 
face resistance is calculated from Eq. 3. Open symbols are 
data obtained with the thinner film 5 = 190nm 
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FIG. 4. Fit of Eq. 2 to a typical curve of Fig. 1. Indicated 
are the separate contributions of three different mechanisms 
of dissipation: 1 is the linear one (which is an exponential 
curve because of the logarithmic abscissa) , 2 is the contribu- 
tion from Eq. 4, and 3 is from Eq. 5. The solid line is the 
sum of these contributions. 

In the nonlinear regime of v{F) the electric field in 
Eq. H must have a nonlinear current dependence. This 
is usually described by thermally activated flux motion 
V{ oc exp{—U{j)/kT) with current dependent pinning en- 
ergies U{j), for a review see . Essentially two different 
dependences are being considered, firstly barriers diverg- 
ing for J — > 0: 



UU) = Uo 



(4) 



where j < jo and < < 1. Vortex glass and collective 
creep models (for the latter ones ^ > 1 is also possible) 
yield this behavior which implies zero dissipation in the 
limit j 0. Secondly, for large currents, one considers 
barriers vanishing at a critical current density jc as 



U{j) = f/o 1 - 



] 



(5) 



which leads to an exponential dissipation for large cur- 
rents (Kim-Anderson model), E oc sinh(j/ji) with ji = 
jckT/Uo- We fit the complete v{F) curves by includ- 
ing both nonlinear mechanisms (Eqs. ^, ^ and a linear 
term additively in the energy balance, see Eq. ^ An 
example is shown in Fig. 4, where the separate contribu- 
tions and their sum are compared with an experimental 
curve. The linear term contributes only at very low am- 
plitudes. In the intermediate range the glass (or creep) 
term (Eq. ^ with — 0.2) is dominant. Only at the 
largest amplitudes the Kim- Anderson term (Eq. |^) be- 
comes relevant. It is this mechanism which leads to the 
observed slow logarithmic increase of the oscillation am- 
plitudes at large drives in Fig. 1 . Towards lower temper- 
atures only the linear term and the glass term contribute 
as no data could be obtained at larger drives because 
the oscillations became unstable. From these fits we can 
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determine the barriers U{j) at various temperatures, see 
Fig. 5. It is a peculiarity of Eq. || that only the constant 
slope —dU/dj — Uo/jc can be obtained reliably from a fit 
to the data because the quantity Uq determines only the 
prefactor which we do not evaluate for it contains other 
unknown quantities, e.g. the density of the trapped flux 
n or the prefactor of the vortex velocity V{. From Fig. 5 
it is evident that at large currents and high tempera- 
tures the barriers of Eq. ^ ultimately drop below those of 
Eq. |4| ijl^. It is interesting to note that all three mech- 
anisms of dissipation are found to apply simultaneously 
and additively. We do not find any indication of a sud- 
den disappearance of one of the mechanisms which might 
be caused by a phase transition in the vortex system. In 
this context it may be important to remember that we 
are dealing with a low vortex density. We do find, how- 
ever, that the critical current jc in Eq. [s] extrapolates to 
zero at 86.5 K (well below Tc = 89.5 K) which might 
indicate an onset of free flux flow where the magnet loses 
its lateral stability. 
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FIG. 5. Dependence of the pinning potential U on the av- 
erage current density in the film obtained from the fits of 
our model to the data: curved lines are Eq. 4, straight lines 
are Eq. 5 at the following temperatures (from top to bottom) : 
78 K, 81 K, 83 K, 84.5 K, 85.5 K. Note that from Eq. 5 we can 
determine only the slopes Uo/jc, in the figure we arbitrarily 
set C/o = 4 • lO^K 

In summary, our method of studying the dynamic 
properties of a levitating magnet yields results on flux 
pinning which are supplementary to those of the more 
standard techniques like measurements of the ac suscepti- 
bility, magnetic relaxation, or current-voltage character- 
istic. Because of the dipolar magnetic field the analysis of 
the results may be more complicated but our method has 
the advantage of being independent of edge effects and 
demagnetization factors. The enormous increase of the 
losses in the YBCO films from a very low level to an ex- 
ponential dissipation is a surprise and is in sharp contrast 



to bulk samples. A more detailed account of our experi- 
ments including an analysis of the results on bulk samples 
(sintered and melt-textured) and of the elastic properties 
of the oscillating magnet will be presented elsewhere. 
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